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Abstract. We describe a new algebraic technique for enumerating self- avoiding walks 
on the rectangular lattice. The computational complexity of enumerating walks of N 
steps is of order 3^/^ times a polynomial in N, and so the approach is greatly superior 
to direct counting techniques. We have enumerated walks of up to 39 steps. As a 
consequence, we are able to accurately estimate the critical point, critical exponent, 
and critical amplitude. 
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1. Introduction 



Over the years, the enumeration of square lattice self-avoiding walks has become 
a benchmark, first for computer performance, and more recently for algorithm design. 
In the early '70s, Sykes et al. (1972) obtained 24 terms of the series by using the chain 
counting theorem. Direct enumeration is probably somewhat faster, but the graphs 
enumerated for the chain counting theorem (figure-eights, theta graphs, dumbells and 
polygons) were useful for other problems in the theory of phase transitions, most 
notably the Ising model. In 1987 Guttmann extended the series by three terms, using 
direct enumeration, and in 1991 Guttmann and Wang using a dimerisation algorithm 
obtained two further terms. Subsequently McDonald et al. (1992) obtained a further 
term, using an extension of dimerisation to trimerisation. Late in 1991, Masand et al. 
(1992) used a large supercomputer, a CM-2 containing 65536 processors, to extend 
the series to 34 terms, running for ~ 100 hours. All these advances came about 
because of improvements in computer technology (in large part), with relatively small 
improvements brought about by algorithm design. (Dimerisation or trimerisation 
saves a factor of around 2 or 3, but does not ameliorate the exponential growth rate 
of computer time). 

The finite-lattice method plus transfer matrices described here allows 35 terms to 
be obtained on a work station (an IBM 6000/530 with 256MB of memory) in less 
time than the 65536 processor CM-2 took to obtain 34 terms. Because of memory 
requirements, it was necessary to move to a larger machine (an IBM 3090/400 with 
500MB memory and 2GB of backing storage) to go from 35 to 39 steps. This computer 
would be capable of extending the series to 43 terms. However that calculation might 
take up to a month of c.p.u. time, and so has not been pursued. This improvement 
however is due to the exponential improvement in algorithm design, rather than 
evolution of computer speed. This is discussed in more detail below. 

The method used is based on the method which Enting and Guttmann have used 
extensively over the past twelve years to enumerate self-avoiding rings on a square 
lattice, but is significantly more complicated due to requirement for a second stage of 
processing. 

In a series of papers we have reported some significant improvements in the 
enumeration of self-avoiding rings on the square lattice extending the known series 
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(Sykes et al., 1972) from 26 steps to 38 steps (Enting, 1980), 46 steps (Enting and 
Guttmann, 1985) and 56 steps (Guttmann and Enting, 1988). The extension from 38 
steps to 56 steps reflects our use of increasingly powerful computing systems and in 
particular the use of increasingly large amounts of physical and virtual memory. Only 
minor changes to our programs have been made, primarily to 'tune' the procedure 
to make efficient use of particular computer architectures. We have also generalised 
our method to enumerate caliper moments of self-avoiding rings on the square lattice 
(Guttmann and Enting, 1988) and also to enumerate self-avoiding rings on the L and 
Manhattan lattices (Enting and Guttmann, 1985) and the honeycomb lattice (Enting 
and Guttmann, 1989). Most recently we have extended the enumeration of triangular 
lattice polygons to 25 steps (as reported by Enting and Guttmann, 1990) and then to 
35 steps (Enting and Guttmann, 1992). 

While our techniques have been highly efficient for enumerating self-avoiding rings 
they are less suitable for enumerating self-avoiding walks. The difficulty is that walks 
can span a larger lattice than rings because they arc not forced to return. A walk of 
L steps can span a distance L while a ring of L steps can only span a distance of up 
to L/2. For walks constrained by a surface our polygon enumeration techniques could 
be generalised in analogy with our calculation of surface susceptibilities for the square 
lattice (Enting and Guttmann, 1980). 

The present paper presents an algebraic technique for enumerating self-avoiding 
walks on a rectangular lattice. The basic quantity that we consider is Cmn-, the number 
of walks from a given origin with n steps in the ±a; directions and m steps in the ±j/ 
directions. We consider segments of walks that double back in the y direction and 
which can therefore be counted efficiently by transfer matrix techniques. The general 
enumeration can be expressed as a combination of such irreducible contributions. We 
further improve the efficiency of the procedure by restricting the range of the index n 
to n <k and reconstruct Cmn for m -|- n < 2A; -|- 1 by using the symmetry relation 

The layout of the remainder of this paper is as follows. Section 2 describes the 
way in which the generating function for self-avoiding walks can be constructed 
from irreducible contributions. Section 3 shows the way in which these irreducible 
contributions can be constructed from generating functions for walks on strips that 
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can be determined by algebraic techniques. Section 4 describes the algorithms 
for determining the requisite generating functions and analyses the computational 
complexity of the procedure. Section 5 describes our analysis of the singularity 
structure of the generating function for self-avoiding walks based on the 39 terms 
that we have obtained. 

2. Generating functions for self-avoiding w^alks 

The generating function for self-avoiding walks on the rectangular lattice is 

oo 

C{U,W)^ J2 CmnUV (2.1) 
m,n=0 

The enumeration of the coeflBcients Cmn is restricted to finite order in m and/or 
n and we generally truncate the double series at m + n < J. The obvious summation 
then gives us the number of walks of up to J steps on the square lattice. 

Our enumeration procedure is based on considering projections of walks onto the 
y axis to produce the type of diagrams shown in Figure 1. We refer to segments of 
walks as irreducible if the projection of that segment onto the y axis has two or more y 
bonds in each position. We will also classify irreducible segments by the number of y 
bonds that they span. At this point we need to refine the terminology and distinguish 
between walks which are directed graphs and chains which are not directed. Our aim 
is to enumerate walks while the transfer matrix techniques generally enumerate chains. 
The decomposition of walks into irreducible components makes use of both chains and 
walks. Figure 2 shows the 5 distinct types of irreducible component that we need 
to consider. 

P{u, w) is the generating function for walks that have no y bonds. Thus 

P{u,w) = l + 2u + 2u'^ + 2u^... = (1 + «)/(!-«) (2.2) 

Q{u,w) is the generating function for chains that are irreducible and for which 
neither end-point lies at an extremal y coordinate. Wc also consider subdividing such 
cases according to m, the number of y bonds spanned by the projection and define 
Qmiu,w) accordingly. The two pertinent results are 
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and 



oo 

Q{u,w) = ^ Qm{u,w) 

m=2 



(2.3) 



QUu,w) = Oiw^^) (2.4) 

Note that from the definition 

Qo = Qi = (2.5) 

R{u, w) is the generating function for irreducible chains with both endpoints having 
the maximal y coordinate. Again we subdivide these chains according to the number 
of y bonds spanned and put 



R{u,w) = J2 Rm{u,w) (2.6) 

m=l 

Rm{u,w) = Oiw^'^) (2.7) 

and arbitrarily define 

Ro = (2.8) 

S{u,w) is the generating function for irreducible chains that have precisely one 
end having the maximal y coordinate. Again 



S{u,w) ^ ^ Sm{u,w) (2.9) 

m=2 

Sm{u,w)=0{w^"'+^) (2.10) 



and from the definition 



So = Si=0 (2.11) 

Finally T{u, w) is the generating function for irreducible chains in which the two 
ends have maximal and minimal y coordinates. We arbitrarily define 

To = (2.12) 
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so that 



oo 



T{u,w)= J2Tm{u,w) 



(2.13) 



m=l 



Tm{u,w) = 0{w^'^) 



(2.14) 



When constructing the generating functions for walks we combine the irreducible 
components represented by P,Q,R,S and T by linking them with single y bonds. 
These each contribute a factor of w to the generating function C{u,w). Two single 
y bonds that are adjacent in the projection are actually connected by a type P 
irreducible component. A chain whose projection is k consecutive single y bonds 
will have a factor of P at each internal point. Note that the walk generating function 
P is required because each distinct direction along the x axis will generate a distinct 
chain when combined with y bonds. If we sum these contributions we can consider 
linking components of types i?, S and T with chains of one or more y bonds connecting 
type P walks. The generating function for such chains is 



In the same way we can regard the overall chain as consisting of end segments P, 
R or S connected by combinations of "reducible" parts with generating function U 
and irreducible parts with generating function T. The generating function for chains 
connecting irreducible end segments is thus 



U{u,w) = w + wPw + wPwPw... 



w/ (1 — wP{u, w)) 



(2.15) 



V{u,w) = U + UTU + UTUTU + ... 



U/{1-TU) 



w/{l-w{T + P)) 



(2.16) 



It is now possible to express the self-avoiding walk generating function as 



C{u, w) = P{u, w) + 2[Q{u, w) + 2R{u, w) + 2S{u, w) + T{u, w)] 

+2V{u, w)[P{u, w) + 2R{u, w) + S{u, w) + T{u, w)f (2.17) 
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The structure of expression (2.17) shows that to obtain an expansion in powers of 
w, it is necessary to obtain Q,R,S, and T to the requisite order. If the irreducible 
generating functions Qm-, R-mi ctnd are known for m < M then C{u,w) will 
be correct to ^y^M+i ^^j^^ g^^g^ incorrect term arising from the absences of Qm+i and 
Rm+i)- The use of the symmetry relations (1.1) will give Cmn for m + n < 4M + 3. 



3. Combining generating functions for chains on strips 

The transfer matrix techniques described in the next section produce generating 
functions for sets of walks confined to strips whose y coordinates are bounded. Subject 
to these constraints, all chains are counted, not merely irreducible components. Thus 
there is a need to relate unrestricted generating functions to the restricted generating 
functions for irreducible components. In our previous enumerations of self-avoiding 
rings only linear combinations of different classes of graph were involved and so the 
restricted generating functions were linear combinations of unrestricted generating 
functions. The present formalism is more complicated because non-linear relations 
are involved. We begin by considering T^{u,w), the generating function for chains 
whose y coordinates span M bonds and which have one end at each y extremum, that 
is, is the generating function for bridges. The sum over M is denoted T*{u,w). 
By considering the appropriate subset of terms from (2.17) we have 

T* {u, w) = T{u, w) + {P{u, w) + T{u, w)fV{u, w). (3.1) 

While this equation is formally correct, it is unsuitable for relating T* to T because 
of the fact that while Tm is of order w^'^, is of order w"^. Thus to obtain T 
correct to would require the calculation of for m < K. This difficulty is 
avoided by introducing an extra variable z whose power corresponds to the width 
of the segment under consideration. We refer to functions including z as "extended 
generating functions". We define the extended generating function for irreducible 
bridges as 

oo 

X{u, w, z) = P(w, w) + ^ z'^Tmiu, w) (3.2a) 

m=l 
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and the extended generating function for all bridges as 

oo 

X*{u, w, z) = P{u, w)+Yl ^'^'^miu, w) (3.26) 

In these terms, V{u, w), the generating function for bridges ending in single bonds, 
generalises to 

V{u, to, z) ~ X{u, w, z) = wz/[l — wzX{u, w, z)] (3-3) 
Relation (3.1) generalises to 

X*{u, w, z) — X{u^ w, z) + X{u^ w, z)'^V{u, w, z) (3-4) 

whence 

X(«, w, z) — X*{u^ w, z)/[l + wzX*{u, w, z)]. (3-5) 

This relation provides the basis of a suitable truncation. The expansion of (3.5) 
to order z^ requires for m = 1 to K and will give for m = 1 to K. The result, 
noted above, that is of order w^"^ provides a useful check on the algebra. 

If we define R^{u, w) as the generating function for chains in a strip of width m 
such that both ends have the maximal y coordinate then 

^ m 

RUu,w) = -{P{u,w) - 1) + J2Rn{u,^) (3-6) 

^ n=l 

This relation can be easily inverted to give individual J?^. These are needed to 
define R{u,w) in (2.17) and also to recover the Sm{u,w) from Y{u,w,z) (eqn 3.8a) 
and the Qm{u, w) from Z(u, w, z) (eqn 3.11) 

We define S^{u,w) as the generating function for chains in a strip of width m 
where one end of the chain has the maximal y coordinate and the other end does 
not have an extremal y coordinate. The general relation between the irreducible and 
unrestricted generating functions is 

S* (u, w) = S(u, w) + [S{u, w) + 2R(u, w)]V(u, w) [P(u, w) + T(u, w)] (3.7) 

We define the extended generating function for irreducible components with one 
or both ends at the maximal y co-ordinate as 
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OO CX) 

Y{u, w,z) = 2Y, z'^Rmiu, w)+J2 ^""Smiu, w) 

m=l m=2 

and the corresponding unrestricted function as 

CO OO 

Y*iu,w,z) = 2 J2 ^"^KnM + z^Sl{u,w), 

m=l m=2 

Eqn. (3.7) generalises to 

Y*{u, w, z) = Y{u, w, z)[l + V{u, w, z)X{u, w, z)] 

or 

Y{u, w, z) = Y*{u, w, z)/[l + V{u, w, z)X{u, w, z)] 

from which the Sm can be recovered once i?„ and T„ are known for n <m. 

FinaUy we consider Q^, the generating function for chains in a strip X 
neither end has an extremal y coordinate. 

We have 

Q*{u, w) = Q{u, w) + [P{u, w) + 2R{u, w)] V(u, w) 

Defining 

OO 

Z{u,W,z) = Y z"^Qm{u,w) 
m=2 

and 

OO 

Z*{u,w,z) J2 z"'Q*{u,w) 

gives 

Z{u, w, z) = Z* (u, w, z) — Y{u, w, tYV{u, w, z) 
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4. Transfer matrix enumeration techniques 



The analysis in the previous sections has reduced the problem of enumerating 
general self-avoiding walks of AK + 3 steps to one of enumerating walks confined 
to strips of width < K subject to various constraints. In order to enumerate walks 
confined to strips we use a transfer matrix technique that generalises the approach that 
we have used in our earlier enumeration of self-avoiding rings. We draw a cross-section 
line (with a kink) across the width of a strip of width K so as to cut + 2 of the bonds 
on which steps of chains can occur. We note that if we specify the set of occupied 
steps then the self-avoidance constraint acts independently to the left and right of 
the cross-section line. However not all combinations of self-avoiding components from 
the left and right of the cross-section line combine to give walks. It is necessary to 
consider the connectivity of the components. This can be done by generalising the 
technique that we used in our earlier work. We assign to each bond intersected by the 
cross-section line an index 



Here "0" denotes an empty bond, "1" denotes a step connected to a (uniquely 
defined) later step, "2" denotes a step connected to a (uniquely defined) earlier step 
and "3" denotes a step not connected to any other steps intersected by the cross-section 
line. 

If we define 



rii = 0, 1, 2, or 3 



1 = 1 toK + 2. 



A{i,j) = {k : k < j and n/. = i} 



(4.1) 



then we require 




for all j 




\A{l,K + 2)\ = \A{2,K + 2)\ 



(4.26) 



as in the enumeration of self-avoiding rings and 



\A{3,K + 2)\ <2. 



(4.2c) 



The numbers of sets of rii subject to these constraints for various K are given in 
Table 1. These numbers give the main limitation on the size of walks that can be 
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obtained because it is necessary to store a partial generating function for the number 
of walks corresponding to each allowed set of n^. These numbers, s^, are larger than 
the corresponding vector sizes, r^, used in the enumeration of self-avoiding walks, but 
only by a factor 7 (A;) which is constrained as 

1 < -f{k) < hk'^ + 5k + 7) 

Thus the increase in the is dominated by a S'^/^ increase as for the r^. 

Self-avoiding chains in strips are developed successively by advancing the cross- 
section line so that one vertex of the lattice passes from the right to the left of the 
line. Except at the beginning of a column, this corresponds to moving the kink down 
one row. This move replaces two bonds (and their associated n^) by two new bonds 
with new n^. The other rii are unchanged except when the addition of the new site 
changes the connectivity of the components. Table 2 shows the various combinations 
of new Ui that can be produced from various combinations of old Ui pairs. Various 
special cases occur. If a link from a free end (i.e. Ui = 3) connects to an existing 
loop segment then the other end of the loop must be reset to type 3. If two loops 
meet then one end must be relabelled. The closing of a single loop implies that a 
ring, disconnected from the walk, has been created and this configuration is ignored. 
The final step in the construction of a chain is when two type 3 bonds meet. When 
this occurs, all other bonds must be empty for a valid chain generating function to be 
added to the running total. If this is not the case it implies that other disconnected 
components are present and the configuration is ignored. 

The iteration is initiated from an empty state (n^ = 0) with generating function 
1. As each bond is added, factors of v or as appropriate are used to multiply 
the old partial generating function and the product is accumulated into the running 
total for the new partial generating function. Each chain could, in principle, be 
generated for a number of different x and y displacements from a given reference 
origin. To ensure uniqueness in the x direction, the chain is required to intersect the 
first column considered. Thus the state with all rii zero is never continued after the 
first column has been built up. This ensures that each chain is counted only once and, 
together with the requirement that the last operation is to join two type 3 bonds, 
also ensures that only one connected component occurs in each graph that is counted. 
Formally, if 4A;+3 series terms are required then the transfer matrix operation must be 
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repeated until 4A; + 3 columns of each strip have been generated. This will ensure that 
all the cancellations involved in going from unrestricted to irreducible contributions 
will be correct. As noted above, the cancellations provide a useful check on the 
implementation of the algebraic formalism. If however the check is not required then 
the results (2.4), (2.7), (2.10) and (2.14) can be assumed to be true. For a strip of 
width M it is sufficient to generate only 4A; + 3 — 2M columns of the strip. 

This requires a large amount of memory to store all the intermediate generating 
functions. If this amount of memory is not available as physical memory, but only 
as virtual (disk based) memory, the performance can be enhanced greatly by being 
careful of the order in which the partial generating functions are processed. This can 
make the entire process close to sequential, and enormously reduce the amount of disk 
access required. 

As mentioned before, when a partial generating function is processed, only the 
two rii coefficients adjacent to the site being added will change, unless this changes 
the connectivity of a loop. However in this case, it will always be a non-zero going 
to another non-zero, so if we make up a "partial signature" out of the rii coefficients 
that are not being directly changed, and just record whether they are zero or not, 
this partial signature will be invariant under the transfer matrix. That is, all "new" 
partial generating functions will have the same "partial signature" . 

We can then process all the partial generating functions with one particular partial 
signature, save the results to disk, and process the next set, until all are processed. 
We do not need to worry about the possibility of having to accumulate a new partial 
generating function to one stored on disk, as if two partial generating functions have 
different partial signatures, they definitely cannot have the same "full signature" {rii 
values) . 

The only problem with this is to make sure that we process all the partial 
generating functions of a particular partial signature first. Fortunately, one can get 
away without having to sort the output file before using it as input. If instead of 
saving to just one output file, two output files are used, it is possible to arrange 
things such that one only has to read from these two files in order, and the data 
will be ordered in exactly the right way for adding the next site. This minimises the 
amount of disk access. To obtain this nice ordering, one looks at whether the new 
value of Hi which will be included in the partial signature for the next phase is zero 
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or not, and accordingly assigns it to one of the two output files. This will arrange the 
partial signatures in binary order, with the most significant bit in the partial signature 
being the most recent bit calculated, and the least significant bit being the oldest bit 
calculated. 

This method has a side benefit — in a multi-processor architecture, one can 
have each processor working independently on separate partial signature groups, with 
very little inter-processor communication. This means that this algorithm is easy 
to parallelise, which will become increasingly important in the future as parallel 
computers are becoming more and more popular, whilst many algorithms are difficult 
to parallelise. 

The problem of ensuring uniqueness in the j/-direction requires inclusion-exclusion 
arguments of the type used in our enumeration of self-avoiding rings. For a strip of 
width M we classify chain generating functions as Gmi^, ±, ±) where -|- denotes an 
allowed location for ends and — denotes a forbidden location. The three arguments 
refer to the top row, the bottom row and the set of internal rows respectively. 
Specifying the end locations permissively, rather than prescribing how many end points 
lie in each set allows us to treat the two ends of the chain independently. We have 

P-1 ^ 

GKi+, -, -) = + J] i?m (4.3a) 

K 

GK{+,+,-) = n + P-l + 2Y,Rm (4.36) 

m=l 

P-1 ^ 

GKi+, -, +) = K^- +Y^(K+l-m){Rm + Q*^ + S*J 

K 

+ Y,{K-m) (T^ + Rm + S*J (4.3c) 

m=l 

K K-1 

GKi-,-,+) = iK-l)^^+J2{K + l-m)Q*^ + 2j2iK-m)iR^ + S*J 

m=l m=l 
K-2 

+ J2iK-l- m)T*^ (4.3d) 

m=l 

These relations can be explicitly inverted to give 
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(4.4a) 



m—l 



Qm = Gm{-, -, +) - Gm-l{+, +) " X/ + Rn + S^) 



(4.46) 



n=l 



P-1 



m — l 



-Qm-E(^n+2*^:+T*) 



(4.4c) 



n=l 



— Gm{+, +, —) — 2Gm(+5 ~) ~) 



(4.4d) 



5. Analysis of series 

This algorithm was implemented for computers in two parts. The first part was 
the transfer matrix portion; that is it computed the G functions. This was the time 
and memory intensive section. It was run up to a width of 8 on an IBM RS6000/520 
with 256MB or RAM, then to width 9 on an IBM 3090 with 500MB RAM and 2GB 
backing store. It required about 200MB and 600MB of memory respectively, and 
required several days in each case. The same machine could have been used to do 
width 10 (43 terms) given several weeks of processor time. This program was written 
in C, and reused memory whenever possible (there was only one bank of memory, 
used for both the about-to-be-processed partial generating functions, and the have- 
just-been-processed partial generating functions). The method described previously 
for using disk storage efficiently was implemented and tested, but not used as whilst 
memory was then no longer a problem, time became a large problem. 

The second part performed all the algebra. Whilst algebra on large (over 15000 
coefficients) polynomials in three variables is slow, it is still a minor problem compared 
to the transfer matrix section, requiring time in only hours and memory in sub- 
megabytes, so efficiency was not as vital. It was implemented in C-I--I-. 

The results for width 9 (39 terms) are given in table 3. 

The method of analysis used is based on first and second order differential 
approximants. It was also used in previous papers [Guttmann (1987), Guttmann and 

15 



Wang (1989) and is described in detail in Guttmann (1989)]. In summary, we construct 
near-diagonal inhomogeneous approximants, with the degree of the inhomogeneous 
polynomial increasing from 1 to 8 in steps of 1. For first order approximants (K=l), 
12 approximants are constructed that utilise a given number of series coefficients, N. 
Rejecting occasional defective approximants, we form the mean of the estimates of 
the critical point and critical exponent for fixed order of the series, N. The error is 
assumed to be two standard deviations. A simple statistical procedure combines the 
estimates for different values of N by weighting them according to the error, with the 
estimate with the smallest error having the greatest weight. As the error tends to 
decrease with the number of terms used in the approximant, this procedure effectively 
weights approximants derived from a larger number of terms more heavily. 

For second order approximants (K=2), we construct 8 distinct approximants for 
each value of N. A summary of the results of this process is shown in Table 4. The 
statistical procedure used to combine the results gives 

Xc = 0.379052 ± 0.000001 7 = 1.3435 ± 0.0003 {K = 1) 

Xc = 0.3790520 ± 0.0000005 7 = 1.3436 ± 0.00015 {K^2) 

These results provide abundant support, if support is still needed, for the value 
7 = 1.34375 obtained by Nienhuis (1982, 1984). To refine the estimate of the critical 
point, linear regression is used. There is a strong correlation between estimates of the 
critical point and critical exponent. This is quantified by linear regression, and in this 
way the biased estimates (biased at 7 = 43/32) are obtained. 

We find 

Xc = 0.3790524 ± 0.0000005 {K = 1) 
Xc = 0.3790525 ± 0.0000005 {K = 2) 

These are in excellent agreement with previous estimates based on the 56 term 
polygon series (Guttmann and Enting, 1988), Xc = 0.37905228 ± 0.00000014 

For the honeycomb lattice, the "connective constant" = 1/xc is known exactly 
(Nienhuis 1982, 1984), and is a/2 + y^, which satisfies a simple quadratic equation in 
x^. A feature of Maple (Version 5) is a clever algorithm for seeking polynomials with 
integer coefficients that have a given root. Attempting to find a quartic polynomial 
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that gave as a root the biased value of Xc quoted above, we found the best solution 
was also a polynomial quadratic in x^. It was 

581a;^ + 7x^ - 13 = 0. 

The root is Xc = 0.37905227.... While we consider it would be fortuitous if this 
were the true value of the critical point, it nevertheless provides a useful mnemonic. 

Another analysis we were able to carry out with this long series was a study of 
amplitudes of the leading term and the correction terms. As previously discussed for 
self-avoiding polygons (Guttmann and Enting 1988), we have found no evidence for 
any non-analytic correction-to-scaling term other than that suggested by Nienhuis, 
with a "correction" exponent of A = 1.5. In the case of the polygon generating 
function this "folds into" the additive analytic term. However, for the SAW series, it 
gives rise to a non-analytic correction term. Furthermore, there is another singularity 
on the negative real axis, at shown by Guttmann and Whittington (1978). 

Thus we expect the generating function for walks to behave like 

C{x) = Ec„x" ~ A{x){l-iix)-^'^^^^[l + B{x){l-iJixf/^...]+D{x){l+iJix)-^'^. (5.1) 

The exponent for the singularity on the negative real axis reflects the fact that this 
term is expected to behave as the energy, and hence to have exponent a — 1, where 
Q! = |. From the above, it follows that the asymptotic form of the coefficients, Cn, 
behaves like: 

Cn ~ //"[ainiV32 ^ «2n-2i/32 + h^n-^'"^^ + {-lyd^n-^'^ + {-iTd^n-'''^] (5.2) 

The five amplitudes, ai, 02, 61, di and d2 come from the leading singularity (giving 
rise to ai and 02), the correction-to-scaling term (giving rise to hi) and the term on the 
negative real axis (giving rise to di and ^2)- A small program written in Mathematica 
was used to fit successive quintuples of coefficients, Cn-4, c^-s, Cn-2, Cn-i and for 
n — 6,7,8,...,39. The results are given in table 5. 

With the possible exception of the sequence {^2}, the sequences for the various 
amplitudes appear to be converging. Various other values for the exponents were also 
tried, including a square-root correction-to-scaling term. In all cases the convergence 
was dramatically worsened by such changes. Indeed, with a square- root correction- 
to-scaling exponent, a number of sequences appeared to diverge rather than converge. 
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However, we have assumed above that the sub-leading term of the singularity on the 
negative real axis is analytic. If we allow this singularity to be a square root singularity, 
so that the last term in eqn. (5.2) above becomes — l)'^(i2n~^ then the results are 
even better converged, as shown in Table 6. 

From these tables we estimate ai ~ 1.1771, a2 ~ 0.554, 6i ~ — 0.19, di ~ —0.19, 
where errors are expected to be confined to the last quoted digit in each case. 
Repeating the above calculations with a critical point shifted by twice the confidence 
limit quoted does not change these amplitude estimates. 

This then completes our numerical study of the generating function for self avoiding 
walks. 
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Tables and table captions 



Table 1. The sizes of vectors required by the transfer matrix formahsm. For ring 
enumeration rk components are required. For walk enumeration Sk components are 
required. 



btrip width, k 


111 1 

#bonds 


rk 


Sk 




1 


1 


2 




2 


2 


5 


1 


3 


4 


13 


2 


4 


9 


37 


3 


5 


21 


106 


4 


6 


51 


312 


5 


7 


127 


925 


6 


8 


323 


2767 


7 


9 


835 


8314 


8 


10 


2188 


25073 


9 


11 


5798 


75791 


10 


12 


41835 


229495 
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Table 2. Allowed transformations of bond indices. For the operations, "Build" 
means incorporate the contribution into the new vector, as defined by the new indices. 
"R(a — > 6)" means apply the change a — > 6 to the other end of the chain in order 
to specify the index in the new vector. "Ignore" means perform no operation, as a 
disconnected ring has been generated. "Acc" means accumulate the vector component 
into the chain generating function if all the other are zero, otherwise the operation 
preceding the "OR" is applied 



Old indices New indices Operation 



{nj,nj+i) 



{nj,nj+i) 



0,0 

0,1 

0,2 
0,3 

1,1 
2,2 

1,2 

2,1 
3,3 

3,1 
3,2 



or (1,0) 

or (2,0) 
or (3,0) 



(0,0), (0,3), (3,0) or (1,2) 
(0,1), (1,0) or (0,0) 

(0,2), (2,0) or (0,0) 
(0,3), (3,0) 
(0,0) 
(0,0) 



Build OR Acc 




R(2 ^ 3) 
R(l ^ 3) 



R(2 ^ 1) 



Build 



Ignore 



(0,0) 



, (1,3) 
, (2,3) 



(0,0) 
(0,0) 



Ignore OR Acc 
R(2 ^ 3) 
R(l ^ 3) 
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Table 3. Numbers of self avoiding walks. 



n Cn 



1 

1 4 

2 12 

3 36 

4 100 

5 284 

6 780 

7 2172 

8 5916 

9 16268 

10 44100 

11 120292 

12 324932 

13 881500 

14 2374444 

15 6416596 

16 17245332 

17 46466676 

18 124658732 

19 335116620 

20 897697164 

21 2408806028 

22 6444560484 

23 17266613812 

24 46146397316 

25 123481354908 

26 329712786220 

27 881317491628 

28 2351378582244 

29 6279396229332 

30 16741957935348 

31 44673816630956 

32 119034997913020 

33 317406598267076 

34 845279074648708 

35 2252534077759844 

36 5995740499124412 



37 15968852281708724 

38 42486750758210044 

39 113101676587853932 
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Table 4. Estimates of the critical point (xc) and critical exponent (7) from first 
order (K=l) and second order (K—2) differential approximants. L is the number of 
approximants used. If L is too small (marked with an 'x'), the estimates are not used 
in the subsequent statistical analysis. 



K 


n 


Xc 


error 


7 


error 


L 


1 


19 


.3790473 


— 


-1.3430184 


— 


Ix 


1 


20 


.3790495 


.0000004 


-1.3432502 


.0001073 


2x 


1 


21 


.3790526 


.0000105 


-1.3435231 


.0016130 


4 


1 


22 


.3790469 


.0000133 


-1.3427806 


.0021290 


4 


1 


23 


.3790468 


.0000071 


-1.3427666 


.0011925 


5 


1 


24 


.3790520 


.0000038 


-1.3436373 


.0006432 


7 


1 


25 


.3790525 


.0000052 


-1.3437004 


.0010102 


9 


1 


26 


.3790508 


.0000050 


-1.3433692 


.0010094 


11 


1 


27 


.3790530 


.0000074 


-1.3437617 


.0013817 


12 


1 


28 


3790519 


0000016 


-1 3435724 


0003713 


12 


1 


29 


.3790519 


.0000010 


-1.3435632 


.0002717 


9 


1 


30 


.3790517 


.0000016 


-1.3435091 


.0004340 


9 


1 


31 


.3790518 


.0000011 


-1.3435230 


.0003276 


8 


1 


32 


.3790514 


.0000017 


-1.3434149 


.0005073 


4 


1 


33 


.3790521 


.0000016 


-1.3436098 


.0004159 


9 


1 


34 


.3790525 


.0000028 


-1.3437270 


.0007447 


11 


1 


35 


.3790518 


.0000003 


-1.3435417 


.0000956 


10 


1 


36 


.3790519 


.0000003 


-1.3435722 


.0001233 


10 


1 


37 


.3790517 


.0000011 


-1.3434730 


.0004568 


8 


1 


38 


.3790518 


.0000009 


-1.3435048 


.0003501 


9 


1 


39 


.3790521 


.0000001 


-1.3436392 


.0000515 


2x 


2 


28 


.3790525 




-1.3437307 




Ix 


2 


29 


.3790520 




-1.3435885 




Ix 


2 


30 


.3790518 


.0000006 


-1.3435431 


.0001628 


2x 


2 


31 


.3790518 


.0000005 


-1.3435300 


.0001315 


3x 


2 


32 


.3790513 


.0000016 


-1.3432041 


.0010873 


3x 


2 


33 


.3790515 


.0000007 


-1.3433885 


.0004434 


4 


2 


34 


.3790519 


.0000003 


-1.3435503 


.0000824 


5 


2 


35 


.3790521 


.0000004 


-1.3436142 


.0001314 


5 


2 


36 


.3790519 


.0000006 


-1.3435607 


.0001972 


6 


2 


37 


.3790520 


.0000001 


-1.3435822 


.0000276 


7 


2 


38 


.3790520 


.0000001 


-1.3435845 


.0000264 


6 


2 


39 


.3790521 


.0000002 


-1.3436174 


.0000616 


4 
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Table 5. Sequences of amplitude estimates. Refer equation (5.2) 



n 


d2 


di 


bi 


02 


ai 


29 


0.0639 


-0.1878 


-0.1999 


0.5584 


1.17700 


30 


0.0666 


-0.1879 


-0.2022 


0.5590 


1.17699 


31 


0.0715 


-0.1881 


-0.1980 


0.5579 


1.17700 


32 


0.0738 


-0.1882 


-0.1999 


0.5584 


1.17700 


33 


0.0781 


-0.1883 


-0.1963 


0.5574 


1.17701 


34 


0.0800 


-0.1884 


-0.1979 


0.5578 


1.17700 


35 


0.0838 


-0.1885 


-0.1947 


0.5570 


1.17701 


36 


0.0855 


-0.1885 


-0.1960 


0.5573 


1.17701 


37 


0.0890 


-0.1886 


-0.1932 


0.5566 


1.17701 


38 


0.0904 


-0.1887 


-0.1943 


0.5569 


1.17701 


39 


0.0936 


-0.1888 


-0.1919 


0.5563 


1.17702 



Table 6. Sequences of amplitude estimates, with the exponent associated with d2 
changed from — | to —2. Refer equation (5.2) 



n 


d2 


di 


bi 


02 


ai 


29 


0.0246 


-0.1902 


-0.2004 


0.5585 


1.17699 


30 


0.0252 


-0.1903 


-0.2017 


0.5589 


1.17699 


31 


0.0266 


-0.1906 


-0.1985 


0.5580 


1.17700 


32 


0.0270 


-0.1906 


-0.1994 


0.5583 


1.17700 


33 


0.0281 


-0.1908 


-0.1969 


0.5576 


1.17700 


34 


0.0283 


-0.1909 


-0.1974 


0.5577 


1.17700 


35 


0.0292 


-0.1910 


-0.1952 


0.5571 


1.17701 


36 


0.0293 


-0.1910 


-0.1955 


0.5572 


1.17701 


37 


0.0301 


-0.1912 


-0.1937 


0.5568 


1.17701 


38 


0.0301 


-0.1912 


-0.1939 


0.5568 


1.17701 


39 


0.0308 


-0.1912 


-0.1923 


0.5564 


1.17702 
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Figure captions 

Figure 1. Schematic representation of projections of self-avoiding walks onto the y 
axis. 

Figure 2. The 5 types of irreducible component from which self-avoiding walks are 
constructed. 

Figure 3. The cross-section line defining the set of lattice bonds which specify the 
partial generating functions. 
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